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Abstract 

We derive an expression for functional determinants in thermal spacetimes as a product over 
the corresponding quasinormal modes. As simple applications we give efficient computations of 
scalar determinants in thermal AdS, BTZ black hole and de Sitter spacetimes. We emphasize the 
conceptual utility of our formula for discussing '1/JV' corrections to strongly coupled field theories 
via the holographic correspondence. 



1 Determinants and quasinormal modes 



In a semiclassical quantisation of gravity the partition function can be written schematically as 



Here g± are saddle points of the Euclidean gravitational action Se- We use g* to collectively denote 



determinants of all the operators controlling fluctuations about the background solution g+. As usual, 
the determinants appear on the numerator or denominator depending on whether the fluctuations 
are fermionic or bosonic respectively. This formula is readily generalised to correlators of operators 
by including appropriate sources in the gravitational action. 

Originally developed with a view to elucidating quantum gravitational effects in our universe [TJ, 
the semiclassical approach to quantum gravity has gained a new lease of life through the holographic 
correspondence [3] . This is because the semiclassical limit in gravity corresponds to the large N 
limit of a dual gauge theory. Almost all works have restricted attention to the leading order large 
N result, the exponent in ([Tj. Exceptions to this last statement include [U (U EJ [7] who studied an 
effect due to the determinant term in pure Anti-de Sitter space. In our recent paper [8] we have 
stressed the fact that the semiclassical determinant term in (pQ), a correction, will be important 

in the framework of 'applied holography'. This is the case because there are interesting effects, in 
for instance thermodynamics and transport, which are either entirely absent at leading order or else 
swamped by physical processes that are not those of prime relevance. 

In general it is technically challenging to compute determinants in nontrivial background space- 
times. Even when the Euclidean saddles are homogeneous the calculation can be involved. For 
example, [9] computed the gravitational determinant on S 2 x S 2 to obtain the semiclassical rate of 
nucleation of black holes in de Sitter spacetime. For most cases of interest, often black hole space- 
times, the eigenvalues of the fluctuation operator cannot be determined explicitly and an evaluation 
of the determinant is difficult without resort to WKB or other approximations (although see for 
instance |10j and citations thereof). 

In this paper we present and derive a formula for determinants in black hole backgrounds, indeed 
finite temperature spacetimes more generally, in terms of the quasinormal modes of the spacetime. 
For (complex) bosonic fields the expression appears in section O equation (|40p . and takes the form 



where z* are the quasinormal frequencies, T is the temperature and e is a simple polynomial 





the metric and any other nonzero fields. The det (— V^J term schematically denotes the product of 
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contribution, fixed by a finite number of local terms discussed in section [U Equation ()4ip gives an 
alternative formulation and (|42p the fermionic analog. 

The quasinormal frequencies of a spacetime determine the time evolution of perturbations of a 
particular field about the background. They are the poles of the retarded Green's function for the 
field in question (see e.g. I12j). In the holographic correspondence they are furthermore the 
poles in the retarded Green's function of the operator dual to the bulk field [13j . While the pole 
with smallest imaginary part determines the thermalisation timescale of the field theory [T3] (within 
a linear reponse regime, c.f. [15j), more generally the poles of the retarded Green's function are 
basic physical quantites in the dual field theory. They could be considered the closest notion to 
well-defined excitations in the strongly coupled theory. Furthermore, over the years a vast effort has 
been invested in the computation of quasinormal frequencies in many black hole spacetimes, for a 
recent review see [16]. Much of the insight gained from these studies can now be translated into the 
physics of one loop determinants. In short, we expect formulae like ([2]) to have the following uses: 

• It allows exact computation of the determinant in certain cases. 

• If a single quasinormal mode contains the physics of interest we can isolate its contribution to 
the determinant. 

• Approximate (for instance WKB) expressions for the quasinormal frequencies give a method 
for approximating the determinant. 

• The quasinormal modes are natural physical quantities in the dual strongly coupled field theory 
which does not have weakly interacting quasiparticle excitations. Expressing the dynamics of 
the field theory in terms of these modes may lead to conceptual insights into strongly coupled 
theories. 

Below we will give a couple of examples of the first of these uses, by (re)computing the one loop 
partition function of scalar fields in thermal AdS, BTZ and de Sitter backgrounds. In our paper 
[8] we used the second fact in the above list; the non-analytic behaviour of a specific quasinormal 
frequency (in the T — > limit) as a function of a magnetic field led to quantum oscillations. 

Our objective in this paper is to give a leisurely exposition of ([2]), together with simple examples 
that illustrate the power of analyticity arguments combined with quasinormal mode sums. Although 
our main interest is in black hole spacetimes, we also consider thermal AdS and de Sitter space as 
simple examples that illustrate the approach. We briefly discuss the numerical evaluation of formulae 
like ([2]) in section [5j In section [6] we note that ([2|) can be generalised to conserved quantum numbers 
other than the frequency, for instance the angular momentum quantum number I. In the discussion 
we mention various possible generalisations and physical applications. 
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2 Warmup examples 



Before deriving a formula for general black hole spacetimes, we will consider some simple one dimen- 
sional examples to illustrate the basic idea. 

2.1 Harmonic oscillator 

The simplest case is a harmonic oscillator. The Euclidean path integral at finite temperature T for 
a complex scalar in one dimension — or in other words for the harmonic oscillator with two real 
degrees of freedom — is 

Z(k) = [ V^V4>e-^ /TdT (\ d ^\ 2 + K W 2 ) , (3) 



where periodic boundary conditions are understood: 4>(t) = <fi(r + 1/T). This reduces to the 
evaluation of a functional determinant: 

Z(k) = N(det(-d 2 T + K)y 1 = nY[(ljI + k)- 1 . (4) 

n 

where N is a normalization factor, the product is a product over eigenvalues, and we introduced the 
bosonic thermal frequencies 

ui n = 2imT , (bosonic) . (5) 

Absorbing an infinite K-independent factor JJ ui~ 2 in iV and using a well known infinite product 
formula gives 

Z(k) = N> — L-= . (6) 
sinh 

To fix N' it is sufficient to note that in the small T or large k limit, Z(k) only gets a contribution 
from the ground state energy Eq: Z(k) ~ e~ E °' T . This correctly identifies Eq = ^/k and implies 
N' = 1/4, so finally 

ZU) = 1 —— = Y e -v^(l+m+n 2 )/T (7) 

4sinh 2 ^ ^ 

which is indeed the 2d harmonic oscillator partition function Z = Tie~ H / T . 

Let us now evaluate the same determinant in a slightly more convoluted way, which however will 
readily generalize to the systems of interest to usQ 

Consider the function Z(k) analytically continued to complex values of k. This function has no 
zeros, and poles whenever the operator — d 2 + k has zero modes, that is, if a solution exists to the 



1 This is essentially an adaptation of an argument in Coleman's lecture on instantons [17]. For another application of 
similar ideas to the problem of computing determinants in AdS/CFT, see [6]. For the general mathematical framework 
to which (a more rigorous version of) this type of reasoning belongs, see for example |18] . 
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Euclidean equation of motion 

{-% + k) 0(r) = O, (8) 

which in addition satisfies the periodicity constraint 4>(t + 1/T) = 0(t). A basis of (not necessarily 
periodic) solutions to the equation of motion is, for any value of k, 

<j> ± ( T ) = e ± ^\ (9) 

(Wick rotated back to real time t = —it and for real k > 0, these become the usual oscillatory 
solutions to the real time equation of motion, with frequency The solutions satisfy in addition 
the Euclidean periodicity constraint if and only if ±y/E = iu n , n G Z, or equivalently if and only 
if sinh (±^r) = 0. Thus Z{k) is a meromorphic function with the same zeros and poles (with the 
same multiplicity) as the meromorphi<J^| function 



/(«) = 77717 ' ( 10 ) 

sinh 2 (^) 



Furthermore we knovcl that when k — > oo, Z(k) — > e~ E °/ T = e~^l T while /(«) — > 4e _v/ ^ //T , where 
the branch of \/k is taken with positive real part. Thus *frjx is a holomorphic function without 
zeros or poles which converges to 1 when k — > oo, that is to say, it must be equal to 1 everywhere. 
This reproduces ((7J). 

2.2 Damped harmonic oscillator 

We now add dissipation to the previous example. This generalisation includes effects that are crucial 
in the black hole case. The partition function for the damped complex harmonic oscillator coupled 
to a 'chemical potential' fi is up to normalization given by |19] 



Z = N\\{{uj n + i^f + 2 1 \uj n \ + K) 1 , (11) 

n 

with 7 > and < m for stability. The normalization factor ./V is taken to be independent of the 
parameters (ft, 7, fJ,) and fixed by requiring the correct, known, undamped 7 = limit. The thermal 
frequencies uo n are still given by ([5]). Notice the absolute value signs around the thermal frequencies 
uj n in flTT]). 



2 Despite the square root there is no branch point at k — because the function is invariant under y/Ti — > —\J7i. 
3 We assume this to be known in the present derivation because the analogous piece of information will be known 

in the black hole case: the large mass expansion of the one loop free energy is determined completely by known local 

expressions. In the present context it can be derived directly from the infinite product representation of the determinant 

by noticing that in the large k limit the logarithm of the product is well approximated by an integral over n. 
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The argument we used for the undamped oscillator was based on complex analyticity, and to 
construct a similar reasoning here, we will need to get rid of the explicit absolute value signs. We do 
this by splitting the product into separate parts Z+, Zq and Z_ containing the terms with positive, 
zero and negative ui n respectively, and treating these factors separately. Let us first consider the 
n > part, analytically continued to complex k: 

Z+{k) = N + Yl ((w„ + iiif + 2 7 w n + k)~ X . (12) 

n>0 

This function has poles whenever a solution to the Euclidean equation of motion 

{-{d T - fi) 2 + 2i 1 d T + K) 0(r) = 0, (13) 

satisfies the Euclidean periodicity and positive frequency constraints. A basis of solutions to the 
equation of motion is given by the two functions 



4>(t) = e Z * T , = fi — ij ± y/n ~ 7 2 ~ . (14) 

(Wick rotated back to real time t = —it, for real values of the parameters and for 7 > 0, these give 
damped oscillatory solutions to the original real time equations of motion, with frequency Re z* and 
damping factor Imz*.) The above solution is periodic in Euclidean time and of positive frequency if 
and only if = ico n , n = 1,2,3,..., or equivalently if and only if F (l + Tpf?) = 00. Thus Z+(k) is 
a meromorphic function of K with the same zeros and poles as Ilz*^) T ^1 + '^ff 
The n < part of Z is 



Z_(k) = N- H {{uj n +i^) 2 -2 1 u n + K) l . (15) 

n<0 

This function has poles whenever a solution to the Euclidean e.o.m. with 7 replaced by —7, 

(-(d T - fi) 2 -2i 7 d T + K) 0(t) = 0, (16) 

satisfies the Euclidean periodicity and positive frequency constraints. A basis of solutions to the 
equation of motion is given by the two functional 



0(t) = e~ s * T , z+ = fi + i~/±v / K-~f 2 + 2i7/x . (17) 

(Wick rotated back to real time t = —it, these give onii-damped oscillatory solutions to the original 
real time equations of motion, with frequency Rez* and damping factor Imf*.) The above solution 



4 The bar notation does not mean the are obtained from the z* by complex conjugation, although it does amount 
to that for real values of the parameters k, 7, fi. 
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is periodic in Euclidean time and of positive frequency if and only if z+ = iu n , n = —1, —2, —3, . . ., 
or equivalently if and only if T (l — = oo. Thus Z^{k) is a meromorphic function of k with the 
same zeros and poles as IIz^k) ^ ~ ) • 
Finally the n = part of Z is 

Z (K)=N (-fi 2 + K y 1 . (18) 



This has poles whenever there is a constant solution to (|T3|) or equivalently to (|16j) . that is whenever 
Q 2 = or equivalently Y\ s = 0. 

Putting everything together, we see that Z(k) is a meromorphic function with the same poles 
and zeros as 

/w=n'(^)nr(-f)n^ e»> 

Z*(«) Z*(«) 

From this we can conclude that 

Z(«) = e p ^/( K ) (20) 

where is a entire holomorphic function of k. To determine P(k), we consider the k — > oo limit. 
A direct computation gives 

log /(«) = - 7 ° g ^ K J + log(4vr 2 ) + (21) 

On the other hand, in the k — >• oo limit the logarithm of the product (jlip is well approximated by 
replacing the sum over n by an integral over n. Introducing a UV cutoff \ui n \ < A and fixing a 
(k, 7, /^-independent normalization by requiring the correct 7 = limit, this gives 

~ T~ ~ 

Comparing these results yields 



log Z{k) = - 7 ° g ^ K; + 0(«-V>) + 0(A~ 1 ) . (22) 



P( K ) = -log(4vr 2 )- 27l °^ ) . (23) 

There can be no corrections in inverse powers of k to this expression since P[n) must be an entire 
holomorphic function, i.e. it has a Taylor series expansion in k with infinite radius of convergence. 

Putting everything together, using T(x + 1) = xT(x) and restricting again to real values of the 
parameters (so the z+ are complex conjugates of the z+) gives the final result 



* T II 



4vr 2 T 



2vrT 



2 



(24) 



The residual dependence on the UV cutoff when 7 7^ is related to the fact that while d K log Z is convergent, the 
same is not true for <9 7 log Z. 



6 



The undamped oscillator result is reproduced from this by putting 7 = and using the identity 
xT(ix)T(— ix) = ir/ sinh(-7rx). 

All of these formulae generalize straightforwardly to systems with multiple quasinormal frequen- 
cies z*, by letting the products range over all quasinormal frequencies. 

2.3 Fermionic oscillator 

The partition function of the undamped fermionic oscillator is 

Z F = JviP ViP e~ ti' T dT ^+ m )V> > (25) 

and hence 

Z F (m) = Ndet AP (d T + m) = N (det A p(-^ + m 2 )) 1/2 . (26) 

Here the subscript 'AP' indicates antiperiodic boundary conditions: vp(r + 1/T) = —iJ)(t). Thus, the 
thermal frequencies are now half integral multiples of 2irT: 

u n = 2-kT '(n + |) , (fermionic). (27) 

Apart from this fact the analysis proceeds exactly as in the bosonic case, giving instead of ([7]) 

Z F (m) = 2 cosh y-^j =e 2T + e 2T . (28) 



For multiple oscillators with mass matrix M we get 

Z F = n2cosh(^) , (29) 

where 21* are the fermionic normal frequencies, obtained by solving (d T + M)ip = for tp = e~ z * T . 

Similarly, if we define the damped fermionic oscillator determinant as the square root of the 
corresponding bosonic one but with antiperiodic boundary conditions, we get instead of ([24 



^^{1¥T -r-r 27T 



3 Euclidean thermal geometries 

We now turn to the problem of actual interest to us, computing determinants of Laplace type 
operators on thermal Euclidean spaces obtained by Wick rotating static spacetimes. 
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3.1 Asymptotically AdS black holes 

Consider a static asymptotically AdS spacetime with a horizon. A specific example is Schwarzschild- 
AdS rf+ i: 

t fa2 = _y( r ) t ft2 + _L^d r 2 +r 2 d f^_ l} ( 31 ) 

2 

where V(r) = 1 — + jp- However the arguments which we will give in the following are not 
restricted to this example, and apply to arbitrary static backgrounds. Geometries with horizons are 
the most interesting as computing determinants in these cases will require treating the positive and 
negative frequency modes separately, similarly to the damped oscillator. 

To compute the one loop partition function at temperature T, we Wick rotate and identify: 

t = -ir, t~t + 1/T, (32) 

so the metric becomes Euclidean and periodic in the Euclidean time r. The fact that the spacetime 
is static means that the metric components are independent of r and that there are no space-time 
mixing terms. The horizon of the original black hole becomes the origin of the Wick rotated Euclidean 
space [20j. The time circle contracts to a point there. For a suitable radial coordinate p, the metric 
near the origin is of the standard flat form 

ds 2 « dp 2 + p 2 d9 2 + ds 2 ± , 9 = 2itTt , (33) 

where ds\ is the transversal metric at the horizon. Regularity of this metric at p = fixes the 
periodicity of the angular coordinate r to be (|32h . and hence the black hole temperature to be T. 
For another suitable radial coordinate z, the metric near the AdS boundary z = is of the form 

L 2 

ds 2 rj —j (dz 2 + dr 2 ) + ds\ , (34) 

where L is the AdS radius. 

For concreteness, say that we want to compute the contribution to the one loop partition function 
in such a background from a massive complex scalar (p. This is given by 

Z(A) = f P0e-^*(- v2+m2 )^ oc -. -i 5T , (35) 

y 1 J det (-V 2 + m 2 ) y ' 

defined with the z — > boundary condition 



~ z 



- A 'I + 0{z 2 )) , (36) 



where A is a root of A(A — d) = (mL) 2 . With this asymptotic behavior, the action is positive 

2 



definite for real A and finite for A > which is the unitarity bound of the dual CFT [21]. The 
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fact that the boundary condition (|36|) is specified in terms of A rather than the mass m 2 will shortly 
motivative us to use analyticity arguments in A rather than m 2 . 

We also require regularity at the origin p = 0. In terms of the complex coordinate u = pe~ %e 
with p and 9 as in (|33p . regularity means that <fi should have a regular Taylor series expansion in u 
and u about the origin. For a mode with thermal frequency quantum number n, we therefore have 
the p — >• behavior 

n = n -ind = u> n /2vT -iu>„r jf > q 

P P " ' (37) 

-_„ = p -n e -tnd _ p -w n /2irT e -iu n T j f n < Q , 

where the w n are the bosonic thermal frequencies as in ([5]). 

Paralleling our warmup examples, we will determine Z(A) by analytically continuing to complex 
A and matching poles and zeros. Because of the nonanalytic dependence on n exhibited in (|37p . 
we will be forced to split the partition function into positive and negative frequency parts, as we 
had to do for the damped harmonic oscillator. In fact much of the following analysis will mimic 
the analysis done for the damped harmonic oscillator. Let us start with the positive frequency part 
Z+(A), defined by (|35p restricted to functions with only positive frequency Fourier coefficients. This 
has no zeros, and poles on complex values of A for which a zeromode exists, that is to say, a solution 
to the equation of motiorjf 

-V 2 + & = , 4>* ~ p- lz * /27TT e~^ T (p -> 0) , (38) 

satisfying the boundary condition (|36p at infinity, which in addition is periodic in r. The latter is 
equivalent to requiring^ 2*(A) = iu n , n = 1, 2, 3, . . ., or equivalently T (l + J^fe) = oo. 

When we Wick rotate the solutions to real time t = —ir, the z+ appear as complex frequencies 
of modes with near-horizon behavior 

oc exp (-iz*(x + 1)) , x = . (39) 

The horizon is at x = — oo, hence these modes satisfy so-called ingoing boundary conditions at 
the horizon as well as (f> ~ boundary conditions at infinity. For physical, real values of A, this 
precisely defines the standard quasinormal modes of the AdS black hole [IB]. Thus we arrive at the 
important conclusion that for the physical, real values of A, the z± are the black hole quasinormal 
frequencies. 



6 The equation of motion near p = implies ~ pT I2 */ 2,rT e Z * T From (|37[1 it follows that positive frequency 
solutions are on the p- tz */ 2wT branch, while negative frequency solutions are on the p+^/a^T hr&nch 

7 It is important to keep in mind that we are considering complex A here. For real A, there will generically be no 
zeromodes; the poles of Z(A) are off the real axis. 
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Stability of the background requires that all quasinormal frequencies have negative imaginary 
part; if not, some modes will be exponentially growing with time. 

For the negative frequency part of the determinant we follow a similar reasoning, except that 
now we need to select the other branch of solutions to the equations of motion. The poles in Z_(A) 
are associated to solutions to the equation of motion with p — > behavior ~ p+«z*/2-7rT e -z*r ^ £ Qr 
which = iu n , n = — 1, —2, —3, . . ., or equivalently r (l — ^f^) = oo. 

After Wick rotating to real time, these solutions satisfy outgoing boundary conditions at the 
black hole horizon, and the z+ are the 'anti'-quasinormal frequencies. Despite what the notation 
might suggest, the sets of these anti-quasinormal modes and -frequencies are not necessarily complex 
conjugates of the standard ingoing quasinormal modes. However for real values of the parameters, 
if the operator V 2 is PT invariant (as is often the case; for example it is true for all examples 
we will consider in this paper, as well as for the planar dyonic black hole in four dimensions we 
considered in [8]), these outgoing modes can be obtained from the ingoing modes by acting with the 
PT symmetry, and then in particular the set of ingoing and outgoing quasinormal frequencies are 
each others complex conjugate. 

The n = contribution can be treated similarly. 

Putting everything together 3, we see that Z(A) is a meromorphic function of A with the same 
poles and zeros as the function Y[ z % V z *z* T (^fy) T ("^f) • Therefore 

z B = e p ( A ) TT ±^p* r ( r ( dh) , (40) 

11 4vr 2 T \2ttTJ \2ttT J ' V ' 

-P(A) is a nonsingular holomorphic function of A, which we will later determine by matching the 
A — > oo behavior. The function P (A) is analogous to P{n) in (|2U[) , To make this more precise, we 
have to address regularization and renormalization, as both the determinant (|35p and the product 
(|4U|) are in fact UV divergent. We will postpone this to the next section, where we will consider 
concrete examples of how this works in practice. We will see that -P(A) is in fact a polynomial in A, 
determined by computing a finite number of integrals over spacetime of local curvature invariants. 
To emphasize the polynomial dependence we will write P (A) = Pol(A). 

As we mentioned, in many cases and are complex conjugates. In this case, the above 
expression reduces to ([2|). 

An alternative and sometimes more useful way of writing this formula is 

zW oKA, n @ n ("+sH»-sr. w 

2* n>0 v 7 v 7 



As with the damped harmonic oscillator for the n = part we can work either with the z* or with the z*, and we 
have PL = FL 2* = FL,^ v 7 ^- 
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where the divergent products over n are understood to be regularized, as is the product over the z* 
(more on this in the next section). 

Although we assumed a single complex scalar in a black hole background, exactly the same 
formula holds for different numbers of bosonic fields. For fermions we get, similar to f|30[) : 



Z F = e™^T\ — : — r~~n rT- (42) 

z* \2 "•" 2-kT) l \2 2wTJ 

or, in product form 

*=^nn("^4)(' + 14)' < 43 > 

2* n>0 v 7 v 7 

Here the are again the quasinormal frequencies, obtained by solving the appropriate Dirac equation 
of motion. 

3.2 Thermal AdS 

Geometries without a horizon have less structure; the quasinormal modes become normal modes, and 
our formalism is correspondingly less useful. However, we begin with such simpler examples as they 
will allow a transparent introduction to the computation of the e Pol ( A ) terms later. Furthermore, 
even in these cases our approach appears to have some technical advantages compared to previous 
computations. 

2 

The metric of AdS^+i in global coordinates is of the form (|3ip with V(r) = 1 + The same 
arguments as for the black hole may be applied to this case. The difference is that the thermal circle 
now nowhere collapse to a point, so we are free to (but of course don't have to) treat positive and 
negative frequencies in one go, like we could for the harmonic oscillator in the absence of damping. 
Exactly the same formulae as for the black hole spacetimes continue to hold, although since there is no 
horizon in the Lorentzian geometry, the Lorentzian interpretation as ingoing and and outgoing modes 
is lost. Because the frequencies z* are real now (as there is no dissipation), we can bring the formulae 
to a form similar to those of the undamped harmonic oscillator using xT(ix)T(— ix) = it/ sinh(-7rx). 
In addition, the frequencies have a direct interpretation as one particle energies (with respect to 
global AdS time): \z+\ = E+. Thus, for a complex boson: 

z s=<= poi<a, nrr!H= ePol(a) n-^. m 

z* 2smh^ Zt l-e~— 

which we recognize as a free multi-particle partition function Z = Tie~ H ^ T . For fermions: 

Z F = e Pol ( A ) ft 2 cosh |^| = e Pol ( A ) \{e^ (l + e^) , (45) 

2* 2* 



11 



which again we recognize as Tie H / T . 

The exact normal frequencies of global AdS^+i are well known. For scalars [22j : 



2n + £ + A 

z n ,i,± = ± 1 , n,£ = 0,1,2,... , (46) 

where L is the AdS radius, £ is the angular momentum quantum number, and A the conformal 
dimension, related to the mass by (ml) 2 = A(A — d). The degeneracy Dg of each frequency equals 
the degeneracy of the £th angular momentum eigenvalue on i.e. 

n (d-i) _ (I + d - 1\ _ (I + d - 3\ _ 21 + d - 2 (I + d - 3\ 
< "V d-1 / V d-1 J' d-2 V d-3 J' 10 

For example when d = 3, this gives Z?£ = 2£ + 1. Somewhat degenerate cases are d = 2, for which 
one should take -Do = 1 and Di = 2 for £ > 0, and d = 1, for which L>o = Z?i = 1 and Dp = for 
£ > 1. The complex boson one loop partition function is then 

Z B = e Pol(A) " E «." 2D< ^^ JJfl-e —J . (48) 

The sum in the first exponential is clearly UV divergent. This is the usual boson zero point energy 
divergence and can be absorbed in the UV counterterms in Pol(A). We will return to this in the 
next section, where we will explicitly compute the full AdS3 partition function. 

3.3 de Sitter space 

The above arguments are not restricted to asymptotically AdS spaces. They apply equally well to 
for example asymptotically de Sitter spaces^ Let us consider pure dSd+i to illustrate this. Static 



patches of this spacetime have horizons that are formally similar to black hole horizons 

2 

The de Sitter metric in static coordinates is of the form (|3ip with V(r) = 1 — j^. It is formally 
obtained from the AdS case by substituting L -4 iL. A causal diamond is given by the coordinate 
patch < r < L. The (cosmological) horizon is at r = L. After Wick rotation the causal diamond 
becomes a sphere. The dS temperature is easily read off as T = ^-j;, by computing the radius of the 
r circle required for regularity. 

The Lorentzian quasinormal modes (ft* have ingoing boundary conditions at the horizon. As for 
black hole horizons, the Wick rotation of this statement is simply regularity for positive frequencies. 
Therefore this is again the relevant boundary condition for the purpose of computing determinants. 

The de Sitter quasinormal frequencies for a scalar are 



Zn,e,± = -i 2n + £ + A± ; n j = 0,1,2,... , (49) 



We thank Dionysis Anninos for his encouragement. 
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and Zn&± = (z n ,e,±)* , where the notation parallels that in (|46p . except that A± are now roots of 
A(A — d) = —irnL) 2 , so in particular for m > d/2, the A± become complex. In AdS, the boundary 
conditions fixed the root A, but because now there is no boundary, both A + and A_ appear. 
Inserting these frequencies in PU|) or (|4"T]) and using the degeneracies Dg from (|4"7|) then produces 
the one loop thermal partition function for a complex scalar in dS^+i, or equivalently the inverse 
determinant of the massive Laplacian on the sphere S d+l . We will do this explicitly in a later section. 

4 Local contributions and examples 

So far we have ignored UV divergences and renormalization, and have hidden them in the as yet 
unspecified term Pol(A), which we claimed is a polynomial in A and obtained from the computation 
of a finite number of local terms. We now turn to these issues. To do this, we will study a number of 
concrete examples, which at the same time will serve to test the above expressions, by reproducing 
some previously derived results. 

4.1 Thermal AdS 3 

We start with the thermal AdS partition function (|48p . For simplicity we will focus on AdS3, which 

2 

has a metric of the form (|3ip with V(r) = 1 + jp, dtti = dcf), with <j> ~ <j> -+- 2tt. We will consider 
a real boson instead of a complex one, which means we should take the square root of (|48p . We 
may combine the products over £, n with multiplicity Dp into a single product over k = 2n + 1 with 
multiplicity k + 1: 

OO a oo 

logZ = Pol(A)- J> + 1)^±^ - 5> + l)log(l-e-^) . (50) 

K=0 K=0 

The last sum converges, but the first one is UV divergent — this is the usual zero point energy 
divergence of a free boson. We can regulate this sum for example by a simple energy cutoff. Recall 
that Pol(A) so far has only been specified to be a polynomial in A which contains UV counterterms 
as well as terms necessary to give the correct large A behavior. We may therefore absorb any UV 
divergence which depends polynomially on A into Pol(A), and write: 

log Z = Pol(A) - J2( K + X ) lo § ( X ~ e -i ^r) . (51) 

K 

The second sum agrees precisely with the finite, nonlocal part of the scalar one loop partition function 
obtained for thermal AdS in |25j by explicitly solving the heat equation and applying the method of 
images. 
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The polynomial Pol(A) is fixed by requiring the correct large A behavior. More precisely, since 
the second term in (|5ip vanishes when A — >• oo, we need 



lim (log Z(A) - Pol(A)) = . (52) 

A— >-oo 

To find Pol(A) we therefore need an independent way of determining the large A behavior of log Z( A). 
Because at large A the scalar becomes very massive, this will be determined by purely local expres- 
sions. These can be found straightforwardly from the heat kernel coefficients of the Laplacian. In 
the next part we briefly recall this technology (for a review, see e.g. |26|). 

4.2 Local terms from heat kernel coefficients 

We may write 

log (A/A ) = - jH j (e- xt - e-^) , (53) 

for some fixed reference Ao, as can be seen by differentiating both sides with respect to A. From this 
observation, we get the following formal heat kernel representation of the determinant of an operator 
D + to 2 : 

/■OO JJ. 

log det(L> + to 2 ) = Tr log(£> + to 2 ) = const. - / — Tr e -*P-Hn 3 ) . (54) 

Jo t 

Here and below "const." should be thought of as some fixed reference determinant. 

For any operator D in n = d + 1 dimensions of the form D = —(V 2 + E), with arbitrary gauge 
and metric covariant derivative V and background field E, one has the following t — > asymptotic 
expansion of the trace of the heat kernel: 

00 

Tie' tD ~ (47rf)-t Y, a k( D ) tk/2 - ( 55 ) 

fc=0 

The heat kernel coefficients have universal expressions in terms of local curvature invariants. If we 
ignore boundary terms, then = for odd k and for even k we have dk{D) = j d n x y/gtr dk(D, x), 
with 1 261: 



a (D,x) = 1, (56) 

a 2 (D,x) = \(QE + R), (57) 
6 

a A (D,x) = — (60V 2 E + 60 RE + 180 E 2 + 12V 2 R + 5R 2 (58) 
360 

-2 R pv R» v + 2 R lwp oR! lvp<J + 30 F^F^ ) . (59) 
General expressions are known up to tiio- 
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Prom this asymptotic expansion we can extract the m — > oo behavior of the determinant, which 
only depends on the heat kernel coefficients oq to a n : 



f°° lit 

logdet^ + m 2 ) = const. - (47r)~fVa fc ( J D) / — e~ tm + ©(m" 1 ) . 

t^L Jo * 



(60) 



k=0 

Indeed, the higher order contributions give finite integrals which manifestly vanish when m — > oo. 
Furthermore, all UV divergent terms are contained in this truncated series. 
The determinant can be regularized as 



poo Ji 

log det A (.D + m 2 ) = const.A - / — f(A 2 t) Tr e - t( - D+m ^ , 

Jo t 



(61) 



where f(x) is some function which goes to zero sufficiently fast when x — > and approaches 1 when 
x — y oo. Therefore f(A 2 t) = 1 except infinitesimally close to t = when A — > oo. Popular choices 
are f(x) = G(x — 1) with the Heaviside step function, or f(x) = (1 — e~ x ) k with k > Using 
the step function, for example in dimension n = 3, this then gives for the regularized version of (|60|) 

logdet A (£> + m 2 ) = const.A + ^ 3/2 (~^~ + 2 ™ 2A ~ ^Y~ m *) ' 

+ (4^rW (" 2A + ' a2(jD) + °( m_1 ) • ( 62 ) 

Here ao = J cPx^/g and for a scalar E = 0, so a 2 = J d?x yfg^R. These terms simply renormalize 
the cosmological constant and Newton's constant. In dimensions n = 4, 5 the couplings associated 
to terms in 04 (D) will also be renormalized, and so on. 

4.3 Result for thermal AdS 3 

Let us now use the result of the previous subsection to determine Pol(A) for the thermal AdS3 
example, according to (|52j) . In AdS3 we have R = —6/L 2 . Using logZ = const. — logdet(Z) + m 2 ), 
m 2 = A(A — 2)/L 2 , and absorbing the (A-independent) A 3 term in the constant, (162p boils down to 
the simple expression: 

\3 f A _ i\2, 



log Z (A) _ const. + / ^ (^L _ (A_£iA^ + 0(A -., . 



(63) 



Pquation (j52|) then fixes, in this regularization scheme, Pol(A) to be exactly the polynomial part of 
this expression. In conclusion, we get our final expression for the exact one loop partition function 
for a neutral scalar in thermal AdS3 in this regularization scheme: 

logZ = const. 4- / AVS (^f - + logn (1 - ^T^' ■ («> 
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where 

= ^' T = SIT' < 65 > 
and "const." indicates that the overall normalization of the partition function has not been fixed; in 
other words the above should be thought of as an expression for log Z/Zq for some fixed reference 

Zq. 

The UV finite local term does not depend on the regularization scheme and agrees with the finite 
local term found in [25]. Note that this is nontrivial: we determined this term by matching the large 
A asymptotics, and the analyticity arguments which we used to derive our determinant formula then 
fixes this to be the exact result valid for any value of A. Note that the analyticity argument does 
not hold for to because, as noted below (|36p . the boundary conditions are not analytic in to, but 
rather in A = 1 ± \J (mL) 2 + 1. And indeed the above expression ([64]) has an infinite series of 1/to 
corrections to the leading large to result. 

For scalars with —1 < to 2 < 0, we can impose either A = A + = 1 + \J {mL) 2 + 1 or A = 



A_ = 1 — -J (mL) 2 + 1 boundary conditions [21]. In the CFT dual, going from A_ to A + boundary 
conditions corresponds to turning on a double trace deformation of the UV CFT [27; in which the 
operator dual to 4> has conformal dimension A_, causing the theory to flow to a new IR fixed point 
in which the operator dual to <j> acquires conformal dimension A+. The central charges of the two 
theories are different. In the bulk this translates to different one loop vacuum energy densities V. 
This difference is obtained by subtracting logZ(A_) from logZ(A + ). In the T — > limit: 

reproducing the results of [HE]. Note that this is finite and independent of the regularization. Since 
the both the difference and the cosmological constant are negative, the A + theory has a smaller 
central charge than the A_ theory, as expected. 

4.4 The BTZ black hole 

We now turn to actual black hole spacetimes. The simplest example is the BTZ black hole [28], for 
which an exact computation of the one loop partition function is possible using our formula, which 
can be compared again to |25] , This provides another test of our formula. The one loop free energy 
of a scalar in a BTZ background was first computed in |54j . 

2 

The BTZ black hole metric with zero angular momentum is of the form (I3ip with V(r) = —M+j^, 
M > 0, d£li = dcj), 4> — </> + 27r. (Note that the thermal AdS3 metric is obtained by putting M = — 1.) 
Regularity at the origin r = r + = L\f~M requires 2ttT = L \j^ - 
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The quasinormal frequencies of a neutral massive scalar of mass m 2 = A(A — 1)/L 2 in a non- 
rotating BTZ black hole background are [551 [2SJ US] 

Zp, s> ± = ±y - 2vrTi(A + 2s) , s = 0,1,2,..., p = 0, ±1, ±2, . . . (67) 

and z PtSi ± = (zp, s ,±)*- Here p is the momentum quantum number along the spatial circle and ± 
distinguishes left and right movers. Note that each p therefore appears twice. Equation (|41h (with 
a factor 1/2 because the scalar is real) becomes 



P 



2-kLT 



log Z = -Pol(A) + ^2 [ lo g (n + 2s + A + i P j + log (n + 2s + A 

s>0,n>0,p 

+1 £ [log (2. + A + i^) + log (2s + A - ^)] (68) 

(69) 



s>0,p 



= -Pol(A) + i £(re + l)log((re + A) 2 + (^) 2 ) , 

where we combined the sums over n and s into a sum over re = 0, 1, 2, . . . with multiplicity re + 1, and 
used log (a; + iy) + log(rr — iy) = log(x 2 + y 2 ). Extracting divergent sums while absorbing as before 
all terms polynomial in A into Pol(A), and using the identity J2 p >i^°S (l + f?) = log sm ^ x = 
tcx — log(7rx) + log (l — e~ 2nx ) , we get 



log Z = Pol(A) + log H (1 - «f +A ) {K+1) . (70) 



K>0 

where 



q = e 2lTif , f = 2mLT = --, (71) 

T 

where r was defined before in (|65p . Since again the second term in (|70p goes to zero when A — > oo, 
and BTZ is locally identical to AdS3, the integrand of the local term Pol(A) will be exactly as in 
the thermal AdS3 case (|64l) . 

The product in (I70j) is related by a modular transformation r — > — 1/r to the product in (I64p . 
in agreement with [25] and as expected on general grounds. This provides a nontrivial check of the 
quasinormal mode formula for the partition function. 

Before moving on we make a brief comment about holography at a quantum level. Notice that 
LT —7- in (|7ip implies q — >• 1 causing log^BTZ to diverge as 1/(LT) 2 . The classical bulk action in 
contrast scales oc LT, so at sufficiently small T, the one loop correction becomes larger than classical 
saddle point contributions and we lose control of the loop expansion. Furthermore, the divergent 
scaling of the correction violates extensivity of the dual boundary theory free energy. However, well 
before this happens, the thermal AdS saddle point comes to dominate the path integral — this is the 
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Hawking-Page phase transition [30J . Similar remarks hold for the large LT limit of ^Ads 3 • Quantum 
gravity comes to the rescue of holography, as one might have anticipated. 
The above results can be generalized to arbitrary quotients of AdS^+il^ 

4.5 dS 2 

In even dimensions the heat kernel large m expansion (|60p involves terms logarithmic in m. One 
might worry that this would produce log A terms and that large A matching thus would force the 
inclusion of log A terms in Pol(A), invalidating our analyticity arguments. To illustrate that this is 
not so, and as another check of our general reasoning, we now consider the case of a scalar in two 
dimensional de Sitter space (which Wick rotates to S 2 ) in some detail. We will consider general 
dSd+i in the following section. 

As reviewed in section (j3.3[) the dS 2 massive scalar quasinormal frequencies are 

k + A± 1 



z K ,± = -i L , A± = -±iv, v = y/(mL) 2 - 1/4 , (72) 



with k = 0, 1, 2, . . . and degeneracy 1. This can be obtained from (|49p by combining k = 2n + £ for 
I = 0, 1. Substituting into (jUJ) with T = produces: 

logZ-Poh» = -2 log(K + A±)+ log(t + n + A ± ) (73) 

n,n>0,± k>0,± 

= - J2 (2r + l)log(r + A±) (74) 

±,r>0 

= ^2C'(-1,A ± )-(2A ± -1)C / (0,A ± ). (75) 
± 

Here we zeta function regularized the series; £ is the Hurwitz zeta function defined by analytic 
continuation of C( s , x ) = Yl n >o( x + n ) _S ' an< ^ we denote ('(s,x) = d s ((s,x). The 'Pol' term should 
be understood as a polynomial in A± or equivalently in v. 

We now extract the large mass asymptotics in order to determine Pol(i^). Using the A — > oo 
asymptotics of the Hurwitz zeta function (see e.g. appendix A of [31| ) 

i ^ (-i)fcj? fc r(fc + s-i) 

^ A ^f(sj^— A^+^i ' (76) 

with Bk the Bernoulli numbers, we find that for v — > oo (i.e. m — > oo): 

log Z - Pol(i/) = (log v 2 - 3) v 2 - ^ log v 2 + 0(1/ v) . (77) 



10 Exact results for such quotients of AdS^+i were also obtained in [7], using a different ('holographic') approach of 
which some elements are somewhat reminiscent of elements in our discussion. 
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On the other hand the large m heat kernel expansion (j60|) gives in two dimensions 

logZ = H s / x ^ (/ ^ e_m2t + t / 7 e ~ m2t ) + ° {l,m) - (78) 

After regularizing by introducing a lower cutoff t > e~ 7 A~ 2 with 7 the Euler constant, and dropping 
an m-independent quadratically divergent term, this can be written as 

logZ = 4^/Vs(m 2 log^r-f log ^) +0(l/m) 

= s y, ^ (-^ iog w + - ? iog tlx?) + o(i/ " ) 

= [ same expression as previous line with A replaced by Arg ] 

+ 5F L ^ log ^ - £ log ^) + ^ • (79) 

Here we introduced an arbitrary renormalization scale Arq- Using i?£2 = 2/L 2 and V0I52 = AirL 2 , 
the first line in the last expression can be evaluated as 

(log u 2 - 1) v 2 - 1 log ^ 2 + (-L - i/ 2 ) log(LA RG ) 2 . (80) 

The second line is a logarithmic renormalization of the couplings p and (f) of the classical action 
5*01. = / d 2 x^fg yjp + 4>Rj. We can drop these terms provided we replace the couplings in the 
classical contribution to the full partition function of the theory by their logarithmically running 
counterparts /o(Arg) arm 0(Arq). 

Comparison between (|77|) and (i80j) thus yields: 

Pol(i/) = 2(1- log(LA RG )) i, 2 + ~ log(LA RG ) . (81) 

D 

Notice in particular that all log v terms have canceled, leaving us indeed with a pure polynomial, 
as needed for our analyticity arguments to be valid. This appears to be a nontrivial test of our 
formalism. 

Summarizing, our final result for the exact, renormalized one loop partition function is 

logZ 52 = Poh» + ^2C / (-l,A ± )-(2A±-l)C / (0,A ± ), (82) 

± 

with A± given by (|72p and Pol(i^) by (|8ip . The determinant of the massive Laplacian on the 2-sphere 
was computed e.g. in [32J (excluding the eigenvalue m 2 of the constant mode and suppressing scale 
dependence). The result is given there in terms of a certain infinite product. At the renormalization 
scale A RG = 1/A our resu lt numerically equals theirs. 
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4.6 dS d+1 

The above computation for dS2 can be straightforwardly generalized to dS^+i (i.e. S d+1 ). Instead 
of (|74j) . we have more generally: 

log Z = Pol - Q( r ) lo S ( r + A ±) ' ( 83 ) 

±,r>0 

where the degeneracy Q(r) is a polynomial in r, obtained from (|47|) and some binomial manipulations 

as 

with £>^ +1) as in ([17|> . For d = 1, we have Q(r) = 2r + 1, for d = 2 this is Q(r) = (r + l) 2 , and so 
on. The conformal dimensions A± are as usual related to the mass by 

d 



A ± = -±iv, v = ^m 2 - (d/2) 2 . (85) 
It is easy to see from the definition of the Hurwitz zeta function that for any polynomial Q(r), we 



have ^2^L (^fjy = Q(—x + 5 S ) £(s,x), where 5 S is an operator acting on a function f(s) aa^ 

Ssf(s) = f(s - 1) . (86) 

Using this observation, we zeta function regularize (|83p as 

logZ sd+1 =Pol(v) + ^2Q(-A±+5 s )('(s,A±)\ s=0 . (87) 
± 

This is just a polynomial in v plus a finite number of standard Hurwitz zeta functions. 

To fix the polynomial Pol(z/), we compare the large v asymptotics of our expression (|87p to the 
asymptotics obtained from the heat kernel coefficients. For the sphere these are all computable, see 
for example appendix A of [33] for some explicit tables. By matching asymptotics, we find that jSj 
is in fact the exact result for odd d+1, that is to say, 

Pol 5d+ i(^) = (d+1 odd). 

Interestingly, the large v asymptotic expansion turns out to be purely polynomial in these cases, 
without any \ jv corrections (only with exponentially small corrections). The first few cases are 

TTU 3 , IT r V 5 V 3 . IT ,U 7 V 5 V 3 - . , 

logics ~ , logZ 9 5~ f 1 ), logZ S 7 ~ — ( j + ). 89 

& b 3 > 12 v 5 3' b 90 v 4 3 270 ; v ; 

Dividing these expressions by the sphere volume gives the (negative of the) renormalized one loop 
free energy density. Comparison to analogous computations of the free energy density in AdS (by 
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1 Foi example, J2T=o ( r + x y = (~ x + C( s > x ) = X ' 2 C( S , x ) ~ 2x((s - 1, x) + £(s -2,x). 
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our methods or as tabulated in appendix A of |6j) shows that the results in the last equation are 
basically the naive analytic continuations of their AdS counterparts. 

For d + 1 even, the situation is slightly more complicated. The polynomial part Pol(z^) does not 
vanish and depends on the renormalization scale A RG - For example for d + 1 = 2 we had (|8ip and 
for d + 1 = 4, we get 

Pol^i/) = (~ + ~ log(LA RG )) ^ + + log(LA RG )) v 2 - log(LA RG ) . (90) 

Furthermore, the large v asymptotic expansion of the full log Z does have an infinite series of \jv 
corrections. 

The result (|87p appears to be a simpler expression than any we have been able to find in the 
literature for the determinant of a massive scalar on a d + 1 dimensional sphere (although it is 
not hard to derive this formula from the conventional expressions for determinants). What made 
things simple here is the constraint of holomorphicity in v, and the fact that unlike the eigenvalues, 
the quasinormal frequencies are linear in the quantum numbers, leading to only elementary zeta 
functions. 

5 Numerical evaluation and WKB 

For many spacetimes of physical interest, including the dyonic black holes we considered in [8], 
closed form analytic expressions for the quasinormal frequencies do not exist, and one has to re- 
sort to numerical computations or approximation schemes. Because our formulae expressing the 
one loop determinant as a product over quasinormal modes are all UV divergent, this poses an im- 
mediate practical problem: How can we extract the UV divergences and obtain the correct, finite, 
renormalized determinants? 

In principle, WKB type approximations can be used to control the UV tails of quasinormal 
mode series to arbitrarily high order in inverse powers of suitable quantum numbers. Examples of 
such computations include \i7\ |4"H], l4"9l [50~1 [5T] . Combined with numerical computations of the low 
lying modes which are outside of the WKB regime, this can be used in principle to compute the 
renormalized one loop determinant to arbitrarily high accuracy. 

A procedure along these lines was developed in [52j for the computation of determinants of flat 
space Laplacians with radially symmetric potentials. In their approach, the computation of the UV, 
high angular momentum (I > L) tail of the determinant is reduced to the evaluation of a number 
of integrals of local expressions similar to heat kernel coefficients, yielding an expansion in powers 
of 1/L, in principle to any desired order. The UV divergences are extracted in analytic form and 
canceled by local counterterms in the usual way. The IR, low angular momentum part (I < L) is 
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computed by applying the Gelfand-Yaglom formula |53t [T7] to each individual angular momentum 
mode, which reduces the computation of the determinant of the radial differential operator at fixed i 
to solving a second order linear ODE. Combining the two parts produces a value for the renormalized 
determinant which becomes exact in the limit L — > oo, and highly accurate already at modest values 
of L provided the UV part is computed to sufficiently high order in 1/L. 

A similar procedure should be possible in principle for efficient and accurate numerical evaluation 
and renormalization of expressions such as (|4ip . Presumably the UV part can again be expressed in 
terms of local quantities, while the IR part must be computed numerically mode by mode. The IR 
part will often produce the qualitatively interesting physics; a number of examples are given in the 
discussion section. 



6 Generalization to other quantum numbers 

If there are other conserved quantum numbers such as a momentum, angular momentum or Landau 
level, we can derive a generalization of equation (I41h where this quantum number takes over the 
role of the thermal frequencies. Let's call the quantum number k, and assume it can take values in 
some discrete index set, k E K. Acting on eigenstates of the operator corresponding to this quantum 
number, the differential operator D of which we wish to compute the determinant can be written 
as D{k), where the part of the operator corresponding to the quantum number k gets replaced by 
a c- number determined by kn The operator D(k) can be thought of as acting on scalars living 
on a dimensionally reduced space, with no explicit appearance anymore of the coordinate dual to k 
(in the black hole case considered previously, this coordinate was r). However, the requirement of 
regularity on the original space can lead to particular additional boundary conditions on the reduced 
space. These could have some nonanalytic dependence on k (in the black hole case we had <p oc p' n '). 
Let us assume they are piecewise analytic on subsets K a oi K, a = 1,2, ... ,N (the analog of the 
positive and negative thermal frequencies). 

Now we drop the constraint k G K and look for pairs (A;*, <p a k^), with k± a to be determined 
complex number, solving the equation of motion 

D(h)^=0, (91) 

subject to the analytic continuation of the boundary conditions appropriate to the subset K a . This 
is the analog of the quasinormal mode equation ([38]) , We denote the solution set of all k± solving 
this problem by K+ jCt . 



As a trivial example, if D = —d T — d x and x takes values on a circle of radius R, the i-momentum k takes values 
in iZandD(fc) = -<9? + fc 2 . 



22 



For generic complex k* not in K, cj)^ will not lead to a well defined function when lifted from 
the reduced to the original space, as it violates the quantization condition. Only when k* G K do 
we have a well-defined function on the original space, and in that case (f)^ is an actual zeromode 
of D. For typical real values of the parameters of the operator D (such as the mass), there will be 
no such zeromodes. However after analytic continuation of these parameters to the complex plane, 
zeromodes k+ E K can be expected to occur at complex codimension one loci in parameter space, 
which are zeros of the determinant. Denoting the degeneracy of the quantum number k by dk and 
by matching poles and zeros as before, we are led to the generalization of f|41 j) : 

det J D = e- pol [] H H (k-h) d \ (92) 
a fc*6-K"*,a keK a 

to be understood in a suitably regularized way, and with Pol a polynomial in the appropriate analytic 
parameters, determined by matching e.g. large mass asymptotics. 

As a trivial application, we can rederive the determinant for S n without knowledge of the quasi- 
normal modes. We take D = — V 2 + m 2 and pick as our quantum number the angular momentum 
I. Then D(l) = 1(1 + n — 1) + m 2 and = — A± defined in (f85|) . Furthermore de = is given 
explicitly in flUJ). Thus 

log det D = -Pol + D £ n) lo §(^ + A ±) • ( 93 ) 

±,e>o 

This is identical to (I83jl , so the remainder of the computation is the same as in the remainder of section 
14.61 Written in this way, it is clear that this can also be formally obtained from log det D = ^2 log A 
by factorizing the eigenvalues X(l), with Pol representing a possible factorization anomaly. 

7 Discussion 

In this paper we have advocated an approach to computing determinants in thermal spacetimes 
using sums over quasinormal modes, combined with analyticity in A arguments to fix a 'local' or 
'UV contribution. We will now discuss possible physical applications of this approach followed by a 
discussion of technical directions in which to build upon our results. 

In applying the holographic correspondence to strongly coupled field theories at finite temperature 
one may occasionally be interested in the contribution of only a few quasinormal modes closest to 
the real axis. Such poles can dominate the low energy spectral functions of the theory. Examples of 
situations where this occurs include hydrodynamic poles [34|, 135]. cyclotron resonances [33 EZ], near 
the onset of bosonic instabilities such as superconductivity [38, 39j HO] and in the vicinity of Fermi 
surfaces |41|, [4"2"1 [4"3"1 14"4"] . It was the last of these that initially motivated us to develop the formalism 
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presented in this paper [8], but we expect all of the examples just listed to lead to nontrivial physics 
at the one loop level. 

Prom the perspective of the holographic correspondence, we have expressed the 1/N correction to 
the free energy of a strongly coupled field theory in terms of the poles of retarded Green's functions 
of operators in that theory. This may suggest a reformulation of the finite temperature theory with 
the poles of Green's functions as the 'elementary' quantities. 

We obtained various exact expressions for scalar determinants. Of these the thermal AdS^ and 
BTZ results matched previously obtained expressions while the compact formula (|87p for de Sitter 
space in arbitrary dimensions, in terms of a finite sum of elementary Hurwitz zeta functions, appears 
new (to us). Similar results can be obtained for thermal AdS in arbitrary dimensions. Furthermore, 
determinants in topological black hole backgrounds should be computable in essentially the same 
way as the BTZ case; the relevant quasinormal modes have been found in [45J. A recent list of other 
spacetimes for which the quasinormal modes are known analytically can be found in the introduction 
of [46]. 

We briefly outlined a scheme for efficient numerical evaluation of renormalized determinants as 
products over quasinormal modes in cases where analytic expressions do not exist, by combining UV 
WKB results with IR numerics. We leave implementation of this to future work. 

Important generalisations of the formalism include the extension to stationary, rather than static, 
spacetimes and to more general asymptotic geometries. This should proceed along the general lines 
sketched for arbitrary quantum numbers in section [6j but some additional care may be needed e.g. 
in doing the Wick rotation or making analyticity arguments. Indeed, at a more formal level it would 
be of interest to prove (or disprove) that the various determinants we are computing indeed have the 
strong analyticity properties we have assumed. 
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